An improved method for intracellular DNA (iDNA) recovery from terrestrial environments

Abstract The simultaneous extraction of intracellular DNA (iDNA) and extracellular DNA (eDNA) can help to separate the living in situ community (represented by iDNA) from background DNA that originated both from past communities and from allochthonous sources. As iDNA and eDNA extraction protocols require separating cells from the sample matrix, their DNA yields are generally lower than direct methods that lyse the cells within the sample matrix. We, therefore, tested different buffers with and without adding a detergent mix (DM) in the extraction protocol to improve the recovery of iDNA from surface and subsurface samples that covered a variety of terrestrial environments. The combination of a highly concentrated sodium phosphate buffer plus DM significantly improved iDNA recovery for almost all tested samples. Additionally, the combination of sodium phosphate and EDTA improved iDNA recovery in most of the samples and even allowed the successful extraction of iDNA from extremely low‐biomass iron‐bearing rock samples taken from the deep biosphere. Based on our results, we recommend using a protocol with sodium phosphate in combination with either a DM (NaP 300 mM + DM) or EDTA (NaP + EDTA 300 mM). Furthermore, for studies that rely on the eDNA pool, we recommend using buffers solely based on sodium phosphate because the addition of EDTA or a DM resulted in a decrease in eDNA for most of the tested samples. These improvements can help reduce community bias in environmental studies and contribute to better characterizations of both modern and past ecosystems.


| INTRODUCTION
Recent technological advancements in DNA sequencing have greatly contributed to our understanding of complex microbial communities in a variety of non-extreme and extreme habitats, including for example, forest soil (Cardenas et al., 2015), permafrost (e.g., Yang et al., 2021), deep biosphere (e.g., Moreno-Ulloa et al., 2019;Vuillemin et al., 2018), and desert environments (e.g., Genderjahn et al., 2021;Schulze-Makuch et al., 2021). However, selecting the appropriate method for DNA extraction is an essential and sensitive step in molecular microbiology as DNA quality and quantity influence the success and reliability of every downstream analysis (Williamson et al., 2011), especially in the field of environmental microbiology (Tanase et al., 2015;Walden et al., 2017).
So far, most of the environmental studies based their interpretations and conclusions on the analysis of total DNA, which is the sum of all DNA pools present in one sample. However, total DNA represents, in general terms, two pools: the intracellular DNA (iDNA), which comes from intact living cells, and extracellular DNA (eDNA), which can persist in the environment and iis comprised of DNA that is not contained within an intact cellular membrane (i.e., dead communities or allochthonous sources; Alawi et al., 2014;Carini et al., 2016). This "dead" fraction can represent more than 40% of the total soil DNA pool (Ascher et al., 2009;Carini et al., 2016). Therefore, molecular biological analyses of DNA based on total DNA extractions can be biased toward higher diversity, causing significant misestimation of single taxa-/taxon-relative abundances leading to misinterpretation of specific environmental processes and their importance for the ecosystem (Carini et al., 2016;Lombard et al., 2011). Accordingly, the separation and extraction of iDNA could help to analyze the living fraction of the environmental sample, especially for low biomass environments where isolating RNA or other activity markers (e.g., intracellular ATP) is particularly challenging.
In this context, Alawi et al. (2014) developed a novel DNA extraction method that simultaneously separates iDNA and eDNA from the same sample. As the iDNA is extracted from intact cells, it provides information regarding the living and potentially active community, whereas the simultaneously separated eDNA bears the background of genetic information derived from dead organisms that can be used for studies on relic DNA (Ibáñez de Aldecoa et al., 2017).
The iDNA and eDNA extraction protocol was initially developed for marine sediment samples but has also been successfully applied to different terrestrial environments, including soils from the Atacama Desert (Schulze-Makuch et al., 2018; and different rock types (limestone, quartz-rich shale, and quartz-rich sandstone) from dryland in Namibia (Genderjahn et al., 2021).
It is important to consider that the iDNA and eDNA extraction procedure requires cells to be separated from the substrate before DNA extraction, leading to lower DNA yields in comparison to total DNA extraction methods which lyse the cells within the sample matrix (Williamson et al., 2011). Previous studies have only addressed DNA recovery problems associated with total DNA extraction from environmental samples and modifications were often required to achieve sufficient DNA yields, especially for low-biomass samples (e.g., Barton et al., 2006;Lever et al., 2015). DNA sorption to soil particles has been successfully reduced in some studies using a highly concentrated phosphate buffer during DNA extraction (Direito et al., 2012;Lever et al., 2015). Additionally, adding EDTA can lead to the chelation of DNA-adsorbing metal ions, thereby allowing for DNA extraction from volcanic ash soils (Rai et al., 2010). This approach has also been successfully used in fluorescence in situ hybridization applications by preventing the adsorption of DNA probes to sediment particles and thereby reducing background noise in the analysis (Morono et al., 2020).
On the other hand, to increase cell recovery, different extraction buffer modifications-including the use of different surfactants-have been applied (Kallmeyer et al., 2008;Williamson et al., 2011).
For example, adding a detergent mix (DM) that includes EDTA, pyrophosphate, and Tween 80 in combination with methanol to a sediment sample has proven to increase cell recovery (Kallmeyer et al., 2008), which consequently can increase iDNA recovery.
However, no further systematic studies on terrestrial samples have been conducted to improve the iDNA and eDNA extraction methods.
Our study combined existing approaches, including using a DM and EDTA to increase cell detachment and DNA recovery from the sample matrix during iDNA and eDNA extraction.
The presented modifications were tested on different soil and rock samples that were taken along the Chilean Coastal Cordillera and aimed to improve iDNA recovery for a more precise characterization of the living community. This new method will help assess ecological questions in future environmental studies, especially when dealing with low-biomass samples (data available in Medina & Horstmann et al., 2023).

| Sample material
During fieldwork in 2019 and 2020, soil and rock samples from different depths were taken from four sites along a climatic gradient in the Coastal Cordillera in Chile Oeser et al., 2018) and were transported on an ice pack and stored at 4°C until processing. The northernmost site-Pan de Azúcar (AZ)-is characterized by arid conditions due to its location within the southern extent of the Atacama Desert. Following the climate gradient to the south, Santa Gracia (SG) reflects a semi-arid ecosystem, La Campana (LC) is located in a Mediterranean climate, and the southernmost site-Nahuelbuta (NA)-is characterized by temperate humid conditions (pictures of all sites in Appendix Figure A1). For each site, a surface and a subsurface sample were processed (e.g., for Nahuelbuta: NA-1 = Nahuelbuta surface sample, NA-2 = Nahuelbuta subsurface sample).
The soils were classified according to the Food and Agriculture Organization of the United Nations (FAO, 2006, see Bernhard et al., 2018 as follows: AZ-Regosol; SG-Cambisol; LC-Cambisol; NA-Umbrisol. Additionally, two samples from drill cores recovered in Santa Gracia  were processed to assess the applicability of the iDNA and eDNA protocol to granitic rock samples from deep biosphere environments. An overview of the processed samples is shown in Table 1.

| Intra-and eDNA separation and extraction
The DNA from all samples was extracted in triplicates using a protocol based on the iDNA and eDNA extraction protocol developed by Alawi et al. (2014). All modifications and alternative setups compared with the latter method are described below ( Figure 1). Three grams of sample material were mixed with 6 mL of 120-mM sodium phosphate (NaP) buffer pH 8.0 (400 mL of 120-mM Na 2 HPO 4 + 29.2 mL of 120-mM NaH 2 PO 4 , sterile filtered over a 0.2-µm PES filter and autoclaved) or with alternative buffers (see below) in sterile 50-mL tubes. Samples were placed on ice for 1 min, followed by agitation on a horizontal shaker for 5 min at 150 rpm.
After cooling the samples on ice for 1 min, the shaking step was repeated. The resulting slurry was centrifuged for 10 min at 500g.
The supernatant containing the eDNA and iDNA fractions was collected in a new sterile 50-mL tube and stored on ice. Another 3 mL of buffer was added to the remaining slurry, followed by an additional round of agitation on the shaker, subsequent centrifugation, and supernatant collection. The entire procedure was repeated twice (total: 4 times). About 15 mL of iDNA-and eDNA-containing supernatant was collected. To separate the iDNA (enclosed in the intact cells) from the eDNA, the extract was filtered through a 0.22µm Sterivex™ filter unit (Millipore) using a sterile syringe, thereby leaving the cells attached to the filter membrane while the flowthrough that contained the eDNA was collected for further extraction. Finally, the filter was rinsed with 3 mL buffer to wash off the remaining eDNA.
To extract the iDNA, the filters containing the cells were cracked open under sterile conditions, and the filter membrane was cut off using a sterile scalpel and tweezers and then transferred to a 2.0-mL screw-cap tube along with glass-and zirconia beads. The iDNA was extracted from the filters using CTAB buffer, phenol:chloroform:isoamyl alcohol (25:24:1), and 10% SDS (modified from Nercessian et al., 2005). After centrifugation at 16,000g and 4°C for 10 min, the aqueous phase was transferred to a new tube and mixed with the same volume of chloroform:isoamyl alcohol (24:1), and the centrifugation step was repeated to remove any phenol residues. The aqueous phase was again transferred to a new 2.0-mL tube and mixed with 1.2 mL (2 volumes) of guanidine hydrochloride (GuaHCl, 6 M in x1 TE buffer pH 6.7). The GuaHCl-mix was transferred onto a silicabased spin column in a sterile 2 mL tube (Zymo Research) and centrifuged at 5000g for 1 min at room temperature. The flowthrough was discarded, and the step was repeated. Then, the silica filter was washed twice by adding 500 µL of washing buffer (50% EtOH, 125 mM NaCl, 10 mM tris, 1 mM EDTA, pH 8.0) and centrifuged at 5000g. To remove any remaining washing buffer, the column was dried via centrifugation at 5000g for 2 min and transferred to a new sterile 2.0-mL tube. Finally, the iDNA was eluted in 100 µL of 10 mM Tris buffer (pH 8.0) via centrifugation at 5000g for 2 min.
For the eDNA extraction, the collected flow-through that contained the eDNA was thoroughly mixed with 30 mL (2 volumes) of GuaHCl, transferred to a sterile 50-mL tube that contained a silicabased spin column with an adaptor for 50-mL tubes (Zymo Research), and centrifuged at 500g for 5 min at room temperature. The flowthrough was discarded, and the step was repeated. Afterward, the column was washed twice by adding 500 µL of washing buffer, centrifuged at 500g, and transferred to a sterile 2.0-mL tube. Then, the column was dried via centrifugation at 5000g for 2 min and transferred to a sterile 2.0-mL tube. Finally, the eDNA was eluted as described above.
A negative control was included in all DNA extractions with all methods: one negative control for eDNA extraction, two negative controls for iDNA extraction, one negative control for the respective filter of the negative control during the separation step, and one extra negative control to check all the reagents during iDNA extraction. F I G U R E 1 Scheme of eDNA and iDNA separation/extraction method. (a) General workflow of eDNA and iDNA separation; (b, c) extraction of iDNA (b) and eDNA (c) following the methods described in this work. Sodium phosphate buffer of 120 mM (NaP 120 mM), sodium phosphate buffer of 300 mM (NaP 300 mM), EDTA buffer of 300 mM (EDTA 300 mM), and sodium phosphate-EDTA buffer of 300 mM (NaP+EDTA 300 mM). Detergent mix (DM).

| Alternative buffers
To systematically analyze and improve the method described by Alawi et al. (2014), we used three different alternative buffers, including negative controls for all of the setups. In addition to the 120-mM NaP buffer from the original protocol, we used a 300-mM NaP buffer (400 mL of 300-mM Na 2 HPO 4 + 29.2 mL of 300-mM NaH 2 PO 4 ), 300 mM of EDTA, and a mix of both NaP and EDTA (300 mM). All buffers were adjusted to pH 8.0, passed through a 0.2µm filter, and autoclaved.
Additionally, we tested all buffers in combination with a DM-as

| Quantification of DNA pools
The concentration of both DNA pools was quantified fluorometrically using the Qubit™ dsDNA HS Assay Kit (Thermo Fisher Scientific).
Additionally, quantitative PCR (qPCR) was performed using a CFX 80°C, 3 s of plate reading). All samples were analyzed by running three technical replicates. A standard with a known concentration (2.5 × 10 8 gene copies) of a 16 S rRNA gene PCR fragment of Bacillus subtilis was used to generate a standard curve via serial dilution (10 1 -10 7 gene copies) and to calculate the efficiency (>90% to <110%) using the BioRad CFX software. Corresponding extraction negative controls and qPCR nontemplate controls were included in the analysis. A melting curve analysis was conducted at the end of each run to identify any nonspecific DNA amplification.

| Cell counting
For the cell counting on the test sample (AZ-2), 3 grams of sample material were extracted with the separation step in the iDNA and eDNA extraction protocol (see above 3.2; Figure 1a), omitting the filtration step. The cells from the "iDNA and eDNA extract" were fixed with 10% Paraformaldehyde in a 1:1 ratio and incubated overnight at 4°C. Afterward, 10 mL of the fixed solution was filtered using a 0.2 μm filter (Anodisc™ 25, Whatman™), and the filters were stored at −20°C.
The cells were stained with SYBR Green I according to the protocol of Noble and Fuhrman (1998) with 0.1% p-phenylenediamine as an antifading agent. Cells were counted using epifluorescence microscopy, using a blue filter set (Leica DM1000 Fluorescence Microscope Filter System I3) and covering 200 fields of view (FOV) or until counting 200 cells (Kallmeyer et al., 2008). The cell number was calculated with the following formula: The transference of eDNA through fine mineral particles to the iDNA pool was checked by spiking 1 µg of genomic DNA into 3 g of DNA-free combusted sand (4.5 h at 450°C, washed twice with NaP buffer, centrifuged for 10 min at 500g and the supernatant was removed as described by Alawi et al., 2014) and the procedure described in Section 3.2 was performed. To compare the recovery rate of iDNA and eDNA, the DNA-free combusted sand was also extracted without adding DNA. The extractions were done in triplicates.
Additionally, 500 ng of genomic DNA was spiked into 3 g of the test sample (AZ-2), and then the soil material was extracted using the selected setups: NaP 120 mM, NaP 300 mM + DM, and NaP + EDTA 300 mM. The recovery was compared with the samples extracted without added DNA. All tests were performed in triplicates.

| Data preparation and statistical analysis
The gene copy number was calculated based on the standard curve and divided by the amount of soil/rock samples in grams (gene copy numbers per gram sample). The mean of the three replicates and the standard deviation were calculated for visualization. An improvement factor was calculated for all samples by dividing the yields of the alternative methods by the yields of the original method, with values above 1 indicating an improvement compared with the method developed by Alawi et al. (2014). To visualize improvements in samples from which no DNA could be extracted using the buffer NaP MEDINA CARO ET AL.
| 5 of 24 120 mM (e.g., SG-34, SG-59), gene copy numbers of the corresponding negative control and the detection limit of the Qubit™ dsDNA HS Assay Kit (0.01 ng µL −1 ) were used. Additionally, a spiked/nonspiked ratio was calculated for the spiked experiment by dividing the values (raw qPCR data or gene copy numbers) of the spiked samples versus the non-spiked samples. Values above 1 indicate an increase in DNA recovery due to adding DNA. All plots were generated with R statistical software (R Core Team, 2022) using the ggplot2 package (Wickham, 2016). Statistical analyses were performed using a Friedman test to check significant differences between the modifications and the original method. p Values below 0.05 were considered statistically significant.

| RESULTS
3.1 | DNA recovery with alternative buffers and the addition of DM on the test sample AZ-2 To optimize iDNA recovery during iDNA and eDNA extraction, different modifications in the separation step of the extraction protocol were tested on the Pan de Azúcar subsurface sample (AZ-2).
This low biomass sample (3.0 × 10 3 cell g −1 , Appendix Table A5)  The measured DNA concentrations for the subsurface sample of AZ-2 were mainly below the detection limit of 0.01 ng µL -1 (Appendix Tables A3 and A4), regardless of the buffer setups used.
Only the eDNA extracted with the NaP 120 mM and 300 mM buffers recovered measurable amounts of DNA (0.016 ng µL −1 with NaP 120 mM and 0.020 ng µL −1 ; Appendix Tables A3 and A4). Therefore, only 16 S rRNA gene copy numbers were considered when evaluating the different extraction setups ( Figure 2).
Extracting with the NaP 120 mM buffer, Alawi et al. (2014) led to iDNA yields of around 7.5 × 10 3 gene copies g −1 soil (Appendix Table A4). Generally, the change of extraction buffer and the addition of DM improved iDNA yield for most of the alternative setups. The most significant improvement compared with the NaP 120 mM buffer was achieved using the NaP 300 mM buffer with the DM, which reached 5.4 × 10 4 gene copies g −1 soil-seven times higher than with the NaP 120 mM method ( Figure 2). The NaP + EDTA 300 mM buffer improved iDNA yields by 3 times without the DM (2.2 × 10 5 gene copies g −1 soil) and up to 5 times with the DM (3.8 × 10 5 gene copies g −1 soil; Figure 2). A slight improvement was achieved when adding the DM to the NaP 120 mM buffer increasing gene copy numbers from 7.5 × 10 3 to 1.2 × 10 4 . The use of the NaP 300 mM buffer without DM (no DM) achieved values of 6.9 × 10 3 gene copy numbers, slightly below the yields of the NaP 120 mM buffer. The setups with only EDTA as extraction buffer resulted in a loss of iDNA, decreasing to around 2.3 × 10 3 gene copies g −1 soil without DM and even further decreasing to 1.5 × 10 3 gene copies g −1 soil after adding the DM ( Figure 2, Appendix Table A1).
The eDNA yield from AZ-2 only improved after using the higher concentrated NaP 300 mM buffer reaching values of 1.6 × 10 5 gene copies g −1 soil which is twice as high as the NaP 120 mM with values around 7.1 × 10 4 gene copies g −1 soil. While the NaP 300 mM buffer with DM gave similar eDNA yields to the NaP 120 mM setup of 4.9 × 10 4 gene copies g −1 soil the eDNA yield decreased for all other setups. Extraction setups containing EDTA (EDTA 300 mM and NaP +EDTA 300 mM) decreased eDNA recovery by at least one order of magnitude ( Figure 3, Appendix Table A2).
F I G U R E 2 16S rRNA gene copy numbers of iDNA from the Pan de Azúcar soil subsurface sample (AZ-2) using four different buffers: a sodium phosphate buffer of 120 mM (NaP 120 mM, according to Alawi et al. (2014) in gray), a sodium phosphate buffer of 300 mM (NaP 300 mM), an EDTA buffer of 300 mM (EDTA 300 mM), and a sodium phosphate-EDTA buffer of 300 mM (NaP + EDTA 300 mM). No DM-no detergent mix added, DMdetergent mix added.
Additionally, cell counting was also performed on the AZ-2 sample to verify and compare the iDNA yield with cell numbers (see Appendix Table A5). The results showed x4 higher cell numbers when using the NaP 300 mM buffer plus DM (1.1 × 10 4 cells g −1 ) in comparison to the NaP 120 mM buffer (3.01 × 10 3 cell g −1 ), which is supported by a similar increase in the qPCR data.

| iDNAadditional tests to exclude eDNA cross-contamination
Another important factor to consider is whether there is eDNA cross-contamination, meaning that the iDNA yield can be increased by the transfer of eDNA to the iDNA pool, attaching directly to the filter material or by fine mineral particles that are retained by the filter membrane during the filtration step. Firstly, we analyzed the attachment of eDNA to the filter membrane by adding purified genomic DNA. The results showed that eDNA did not attach to the filter since the signal measured in the iDNA pool was similar to the negative control ( increased up to 7−8 times (except for NaP 300 mM + DM) which indicated that the spiked DNA goes into the eDNA as expected (Table 3, Appendix Table A7). On the other hand, the spiked/ nonspiked ratio for the iDNA pool was around 1, confirming that eDNA is not transferred into the iDNA pool by attaching to the filter material or mineral particles. Only the sample extracted with the NaP 300 mM + DM setup did not show the expected ratios. For the eDNA and iDNA, the ratio was 0.55 and 0.22, respectively, since we found that the spiked surprisingly sample showed a lower DNA recovery in comparison with the non-spiked sample (Table 3).
3.3 | Evaluation of the improved extraction methods on soil and drill core samples

| NaP 300 mM + DM extraction setup
To verify the NaP 300 mM buffer + DM setup, we tested it on a variety of environmental samples, including Nahuelbuta (NA-1 and NA-2), La Campana (LC-2), Santa Gracia (SG-1, SG-2, SG-34, and SG-59), and Pan de Azúcar (AZ-1). DNA concentration and gene copy F I G U R E 3 16S rRNA gene copy numbers of eDNA from the Pan de Azúcar soil subsurface sample (AZ-2) using four different buffers: a sodium phosphate buffer of 120 mM (NaP 120 mM, according to Alawi et al. (2014) in gray), a sodium phosphate buffer of 300 mM (NaP 300 mM), an EDTA buffer of 300 mM (EDTA 300 mM), and a sodium phosphate-EDTA buffer of 300 mM (NaP+EDTA 300 mM). No DM, no detergent mix added; DM, detergent mix added.
T A B L E 2 Spiked experiment with the Buffer NaP 120 mM and sand. eDNA and iDNA yields, with raw qPCR data, of nonspiked and 1 μg genomic DNA spiked with the Buffer NaP 120 mM and free-DNA combusted sand.  Except for the Santa Gracia subsurface sample (SG-2), all samples showed an increase in iDNA recovery reflected by an improvement of either DNA concentration or gene copy number using the NaP 300 mM + DM. This setup was statistically superior (p < 0.05) to the NaP 120 according to the Friedman test (Appendix Figure A3).
Additionally, the deep biosphere sample SG-34 yielded amplifiable DNA of 1.4 × 10 3 gene copies g -1 rock after extraction with NaP 300 mM + DM (Appendix Table A3 and Figure A2). Small improvements with this buffer were observed for both Nahuelbuta samples Compared with the iDNA, eDNA recovery from the samples was generally less improved by the alternative buffers (Figure 4b) as is also shown using the Friedman test where no significant differences were found with the different setups (Appendix Figure A4). However, the NaP 300 mM buffer plus DM (NaP 300 mM + DM) achieved an Note: eDNA and iDNA yields, yields, in 16S gene copy number g-1, of the nonspiked and 500 ng genomic DNA spiked with the test sample AZ2 using "NaP 120 mM". "NaP 300 mM + DM" and "NaP+EDTA 300 mM" setups.
F I G U R E 4 Improvement factors with the 300 mM sodium phosphate buffer with detergent mix (NaP 300 mM + DM) setup. Improvement factors of concentrations and 16S rRNA gene copy numbers for iDNA (a) and eDNA (b) for the different samples: Pan de Azúcar (AZ-1, AZ-2), Santa Gracia (SG-1. SG-2), the deep biosphere (SG-34, SG-59), La Campana (LC-1, LC-2) and Nahuelbuta (NA-1, NA-2). Values > 1 on the x-and y-axes indicate an improvement in gene copy numbers or DNA concentrations compared with the original NaP 120-mM method by Alawi et al. (2014). Samples within the red area showed a decrease in DNA recovery, while those within the blue area showed an increase in both DNA concentrations and gene copy numbers. Scales are logarithmic.
For the Nahuelbuta sample, eDNA increased by a factor of 6,250 ( Figure 4b), or from 4.8 × 10 6 to 3.0 × 10 10 gene copies g -1 soil (Appendix Table A4 and Figure A2). However, this increase is not consistent with the measured DNA concentrations, which dropped from 53 to 0.43 ng µL -1 (Appendix Table A4). For the Pan de Azúcar surface sample (AZ-1), the use of the NaP 300 mM + DM setup resulted in a 50-fold higher eDNA yield (Figure 4b). Moreover, for samples SG-34 and SG-59, which showed no measurable eDNA when using the NaP 120 mM buffer, gene copy numbers of more than 1.5 × 10 3 g -1 rock were achieved using the NaP 300 mM + DM extraction setup (Appendix Table A4 and Figure A2). This corresponds to an improvement of at least 20-fold (Figure 3b), considering the detection limit (i.e., 396 gene copies g -1 soil) of the qPCR run for the SG-34 and SG-59 samples, which had been extracted with the NaP 120 mM buffer.

| NaP+EDTA 300 mM extraction setup
Additionally, we tested an EDTA-based buffer, the NaP+EDTA 300 mM extraction setup without DM, on all samples ( Figure 5).
Using this setup, the iDNA of the test sample AZ-2 increased from 7.5 × 10 3 to 2.2 × 10 4 gene copies g -1 soil (3 times higher, Figure 5a, Appendix Table A3 and Figure A2) which was also shown by the increase in cell numbers from 2.97 × 10 3 to 2.82 × 10 4 (around nine times higher, Appendix Table A5).
Using the NaP + EDTA 300 mM extraction setup, the Pan de Azúcar surface sample (AZ-1) showed a 4-fold iDNA improvement compared with the NaP 120-mM method and increased from 5.8 × 10 3 to 2.1 × 10 4 gene copies g -1 soil (Figure 5a). Moreover, the Nahuelbuta subsurface (NA-2) sample doubled iDNA yields ( Figure 5a) from 1.8 × 10 7 to 4.7 × 10 7 gene copies g -1 soil. This difference became even more apparent in the measured iDNA concentration, which rose more than 5 times from 1.16 to 6.10 ng µ L -1 (Figure 5b, Appendix Table A3, and Figure A2). An important improvement with this extraction setup was achieved for the two deep biosphere samples from Santa Gracia (SG-34 and SG-59), for which no DNA could be extracted using NaP 120 mM. These samples yielded enough iDNA to be amplified and quantified by qPCR with values of around 2 × 10 3 gene copies g -1 rock (Appendix Table A3 and Figure A2). This shows that the NaP + EDTA 300-mM buffer led to the best extraction results for the deep biosphere samples analyzed in this study.
For the NaP+EDTA 300-mM setup, only the eDNA from the Pan de Azúcar surface sample (AZ-1) and the Santa Gracia surface sample (SG-1) showed improvements. For AZ-1 eDNA, recovery was 4-fold higher (Figure 5b), increasing from 9.1 × 10 3 to 3.8 × 10 4 gene copies g -1 soil (Appendix Table A4 and Figure A2), whereas for SG-1, a 3-fold improvement was observed (Figure 5b), with gene copy numbers of up to 1.3 × 10 7 g -1 (Appendix Table A4 and Figure A2). Further, the separate extraction of the iDNA and eDNA pool from the same sample is a unique approach that has thus far only been F I G U R E 5 Improvement factors with the 300 mM sodium phosphate buffer plus 300 mM EDTA (NaP + EDTA 300 mM) setup. Improvement factors of concentrations and 16S rRNA gene copy numbers for iDNA (a) and eDNA (b) for the different samples: Pan de Azúcar (AZ-1, AZ-2), Santa Gracia (SG-1. SG-2), the deep biosphere (SG-34, SG-59), La Campana (LC-1, LC-2) and Nahuelbuta (NA-1, NA-2). Values > 1 on the x-and y-axes indicate an improvement either for gene copy numbers or for DNA concentrations compared with the original NaP 120-mM method by Alawi et al. (2014). Samples within the red area showed a decrease in DNA recovery, while those within the blue area showed an increase in both DNA concentrations and gene copy numbers. Scales are logarithmic.
applied in a few studies in extreme terrestrial ecosystems (e.g., Genderjahn et al., 2021;Schulze-Makuch et al., 2018. In this approach, the most critical processes that influence its success are the detachment of cells and the desorption of DNA from the sample matrix. In this context, cell detachment from the sample matrix requires that fundamental parameters, such as physical and chemical disper- shaker, using four washing steps that could recover over 90% of the eDNA as described by Alawi et al. (2014).
In terms of chemical dispersion, however, no clear consensus has yet been reached, and various specific extraction buffers have been described (Lever et al., 2015;Williamson et al., 2011). Support has been found for using phosphate buffers (Alawi et al., 2014;Corinaldesi et al., 2005;Lever et al., 2015;Ogram et al., 1987), and particularly for highly concentrated phosphate buffers (Direito et al., 2012;Lever et al., 2015). In addition, using EDTA has successfully prevented DNA adsorption (Rai et al., 2010

| Cell detachment: iDNA yields
Various mechanisms keep cells attached to surfaces, such as direct adhesion (Corinaldesi et al., 2005) and biofilm formation, which are involved in the colonization of mineral surfaces (Costa et al., 2018;Fletcher & Floodgate, 1973). Therefore, a major challenge in successfully recovering iDNA involves the detachment of cells from the sample matrix while maintaining cell integrity. Adding a DM in the washing steps successively increased iDNA recovery. This finding is consistent with observations in which the DM contributed to increased cell detachment for quantitative analysis, such as cell counting (Kallmeyer et al., 2008). Notably, our results indicate that the use of a DM has a variable effect on iDNA recovery depending on whether it is used in combination with a phosphate buffer (NaP 120 mM, NaP 300 mM), an EDTA buffer (EDTA 300 mM), or a mixture of both (NaP+EDTA 300 mM). For most of our results, we observed an increase in iDNA recovery when the DM was added. The most remarkable improvements occurred when using the detergent mix in combination with sodium phosphate (NaP 120 mM, NaP 300 mM) or when using sodium phosphate mixed with EDTA 300 mM (NaP + EDTA 300 mM). This indicates that phosphate constitutes an important part of the functionality of the DM.
Detergents are commonly used in total DNA extractions to lyse cells for a higher DNA recovery (Lever et al., 2015). However, cell lysis should be avoided for iDNA and eDNA extraction. For example, the addition of low amounts of SDS has been reported to cause cell lysis for several methanogenic archaea (Wagner et al., 2013;Wagner, 2020;Wu & Lai, 2011), which would move DNA from the intracellular to the eDNA pool and create bias (Alawi et al., 2014).
However, our data revealed no systematic indication that the DM lyses cells. Therefore, we assume that low amounts of Tween 80 and sodium pyrophosphate in the DM do not lead to cell lysis before separating cells from eDNA in the extraction procedure.
Although phosphate buffers are widely used in DNA extraction protocols (e.g. Alawi et al., 2014;Corinaldesi et al., 2005;Lever et al., 2015;Ogram et al., 1987), no study has yet reported on the relationship between phosphate concentrations and cell detachment processes. The implementation on the test sample from Pan de Azúcar (AZ-2) revealed that the sole use of NaP 300 mM (without DM) did not improve iDNA recovery. This means that high phosphate concentrations in the extraction buffer do not increase cell detachment in our study.
Pure EDTA 300 mM buffer without DM showed a strong decrease of iDNA recovery for the low-biomass sample from Pan de Azúcar (AZ-2, Figure 2). This indicates that EDTA 300 mM alone is not appropriate to recover cells from inorganic, low-biomass samples (e.g., AZ-2) that contain, for instance, notable amounts of metal ions . However, the extraction of iDNA improved substantially after applying the combination of the phosphate buffer and EDTA 300 mM. Our result suggests that EDTA can have a beneficial effect on cell detachment when it is used together with a phosphate buffer. Little knowledge exists thus far on the effect of EDTA on cell detachment processes, and the role of EDTA in this process cannot yet be fully assessed. For example, certain organisms can attach to the mineral surfaces of iron oxides (Neal et al., 2005).
The dissolution of iron oxides by EDTA (Nowack and Sigg, 1997) could lead to the release of cells and could, therefore, improve iDNA recovery. Additionally, one of the main components of biofilms is represented by extracellular polymeric substances (EPSs). These EPSs contain divalent cations (e.g., Mg +2 and Ca +2 ), which are essential to the integrity of the biofilm (Nielsen & Jahn, 1999). EDTA can detach cells from the matrix because it destabilizes the EPS by chelating divalent cations (Nielsen & Jahn, 1999), consequently increasing iDNA recovery.
Our alternative extraction setups successfully improved cell detachment as indicated by a correlation between increasing iDNA yields and cell counts for the modified methods on AZ-2. While it is often difficult to compare DNA yields to direct cell counts, both measurements show similar improvement ratios with the updated methods. Furthermore, the spiked experiments showed no traces of added eDNA in the iDNA pool. This confirms that there is no crosscontamination of eDNA on the filter membrane during the separation step. Particularly, in the spiked experiment using AZ-2 with the NaP 300 mM + DM setup, the eDNA yield did not show the expected increase since the ratio was below < 1, we attributed this due to the effect of the DM which was involved in the reduction of eDNA yields.
However, iDNA yields also decreased after adding genomic DNA by 78%. This reduction might be explained because of sample variability between the two extractions since 16 S rRNA results tend to vary more between different technical replicates, especially for low biomass samples like AZ-2 (Claassen-Weitz et al., 2020).
Summarizing the results of the iDNA yields of AZ-2, the combination of highly concentrated phosphate buffer in combination with either DM or EDTA could improve cell detachment and, thus, increase iDNA yields.

| DNA desorption: eDNA yields
In soils and sediments, eDNA is often adsorbed to particles by binding to -OH groups of for example iron oxides like ferrihydrite, negatively charged surfaces such as clay minerals, or organic matter (Saeki & Kunito, 2010). It is more challenging to extract when the eDNA is attached to mineral and soil particles. Therefore, maximizing the desorption of DNA molecules from mineral surfaces and organic compounds is critical to extract eDNA efficiently. The NaP 300 mM buffer without DM (NaP 300 mM) was the best extraction buffer for the AZ-2 test sample. Large amounts of phosphate in extraction buffers are known to improve total DNA yields because phosphate competes with the negatively charged phosphate backbone of the DNA molecule for binding sites at charged mineral surfaces (Ogram et al., 1987). However, since the DNA adsorption capacity of the sample material is limited, an excess of phosphate due to highly concentrated buffers could be carried over and then compete with the -OH groups of the DNA for the silica membrane in the elution step, reducing the DNA yield (Lever et al., 2015;Lloyd et al., 2010). In addition, an excess of phosphate has also been suggested to lead to the co-extraction of organic compounds (e.g., humic acids), which could also affect DNA recovery and inhibit DNA-based analyses, such as PCR (Lever et al., 2015). Overall, phosphate concentration is an important aspect to consider in terms of the adsorption capacity of the sample, which can only be determined empirically (Lever et al., 2015). However, our data shows no indication of decreased eDNA recovery when using a higher concentrated sodium phosphate buffer exclusively (without DM).
High buffer concentrations often require sample incubation at high temperatures of over 60°C to prevent the precipitation of salts during DNA extraction. Incubations at high temperatures and any other major temperature fluctuations must be avoided in the iDNA and eDNA extraction protocols to ensure cell integrity during the washing step. Therefore, we kept our samples on ice during incubation. Although minor visible salt precipitation occasionally occurred while extracting with NaP 300 mM, iDNA and eDNA recovery were not negatively affected when using this setup, making 300 mM of sodium phosphate in the buffer suitable for extracting iDNA and eDNA. This aligns with Lever et al. (2015), who suggested using 100-1,000 μmol of PO 4 g -1 soil for organic-rich sediment samples and even higher concentrations for organic-poor clay sediments due to their high sorption capacity. These values are comparable to the concentrations used in the present work, which were~320 μmol of PO 4 g -1 soil in the NaP 120 mM buffer and~800 μmol of PO 4 g -1 soil in the NaP 300 mM buffer.
The low eDNA yields when applying the EDTA 300 mM buffer to AZ-2 indicate a negative effect of EDTA 300 mM on eDNA recovery.
However, a concentration of 300 mM of EDTA has previously been used to improve total DNA desorption from volcanic ash soils that contained amorphous aluminum (Rai et al., 2010). The negative effect we observed with EDTA might be explained by the specific composition of our samples or by the possible co-extraction of inhibitors. Consequently, when EDTA and sodium phosphate were used in the NaP+EDTA 300-mM buffer, we observed a substantial decrease in eDNA recovery compared to NaP 300 mM alone. Thus, phosphate buffers are the most efficient agents for desorbing eDNA from the sample matrix.
Moreover, the precipitation of salts increased significantly after adding the DM to the sample buffer mixture because the DM contained 100 mM of sodium pyrophosphate. Precipitation during the washing step could explain the decrease of iDNA and eDNA yields after adding DM in some setups. Notably, such precipitation and its negative effect on DNA yields occur only occasionally and are therefore likely dependent on the specific sample composition (Lever et al., 2015).

| Evaluation and method validation with different soil samples
Based on the tests on the low biomass sample from the subsurface of Pan de Azúcar (AZ-2) we identified an alternative protocol that significantly improved iDNA recovery as reflected in the increase of gene copy numbers g -1 and cell g -1 soil while still providing sufficient eDNA yield for further downstream analyses. We tested the most promising NaP 300 mM buffer plus DM on various terrestrial samples that represented environments with different climate conditions and that also covered both surface and subsurface soils. With this setup, it was possible to increase iDNA for most of the analyzed samples (Friedman test p < 0.05, Appendix Figure A3).
Extreme environments often exhibit rare geochemical compositions, including high concentrations of metals or salts. Since the use of EDTA and phosphate buffer was successfully applied during DNA extractions from particular samples containing amorphous aluminum (Rai et al., 2010) and improved the iDNA recovery on the test sample AZ-2, we also tested the NaP + EDTA 300 mM setup without DM on different terrestrial samples. This setup was applied to the same set of environmental samples as the NaP 300 mM buffer with the DM.
The NaP+EDTA 300 mM setup was not statistically significant in comparison to the NaP 120 mM setup (Friedman test p = 0.342, Appendix Figure  weathering-induced iron oxides, adding EDTA 300 mM as a chelating agent can be beneficial for iDNA extraction. On the other hand, for most of the samples, the best eDNA results were archived using pure phosphate buffers without additives (NaP 120 or 300 mM) because the addition of EDTA 300 mM and DM resulted in decreased eDNA concentrations. Nevertheless, the best eDNA result for deep subsurface samples SG-34 and SG-59 was achieved using the NaP 300-mM + DM setup, again indicating the DM variable effect.
Our study aimed to improve the iDNA recovery using the iDNA Access funding enabled and organized by Projekt DEAL.
Buffer Detergent mix addition Replicate 16 S rRNA gene copy number (Gene copies g-1) NaP 120  T A B L E A2 eDNA yields from Pan de Azúcar subsurface sample (AZ-2) using sodium phosphate buffer 120 mM (NaP 120 mM), sodium phosphate buffer 300 mM (NaP 300 mM), EDTA buffer 300 mM (EDTA 300 mM) or sodium phosphate-EDTA buffer 300 mM (NaP+EDTA 300 mM) extraction setups. T A B L E A3 iDNA yields from all samples using the NaP 120 mM, NaP 300 mM +DM, and NaP + EDTA 300 mM setups. T A B L E A4 eDNA yields from all samples using the NaP 120 mM, NaP 300 mM +DM, and NaP+EDTA 300 mM setups. Note: Missing values due to measuring errors or outliers are indicated as NA. As the improvement factors are based on the relative improvement compared to the NaP 120 mM, improvement factors for this setup are indicated as NA. For calculation and visualization of improvements for samples that were not measurable with the NaP 120 mM method, the detection limit of the Qubit fluorometer or the respective qPCR run was used. These values are marked in italicin this table. Detergent-mix (DM).
T A B L E A5 Cell counting of the test sample AZ-2 using the "NaP 120 mM". "NaP 300 mM + DM" and "NaP+EDTA 300 mM" setups. Note: eDNA and iDNA yields, yields, in 16 S gene copy number g-1, of the non-spiked and 500 ng genomic DNA spiked with the test sample AZ-2 using "NaP 120 mM". "NaP 300 mM + DM" and "NaP+EDTA 300 mM" setups.
F I G U R E A1 Picture of the location of the samples used in this study. (a) The northernmost site-Pan de Azúcar (AZ)-is characterized by arid conditions due to its location within the southern extent of the Atacama Desert. (b) Following the climate gradient to the south, Santa Gracia (SG) reflects a semi-arid ecosystem, (c) La Campana (LC) is located in a Mediterranean climate, and (d) the southernmost site-Nahuelbuta (NA)-is characterized by temperate humid conditions. F I G U R E A2 16S rRNA gene copy numbers using four different buffers: a sodium phosphate buffer of 120 mM (NaP 120 mM, according to Alawi et al., 2014), a sodium phosphate buffer of 300 mM (NaP 300 mM), and a sodium phosphate-EDTA buffer of 300 mM (NaP + EDTA 300 mM). 16 S rRNA gene copy numbers for iDNA and eDNA for the different samples: Pan de Azúcar (AZ-1 and AZ-2), Santa Gracia (SG-1 and SG-2), the deep biosphere (SG-34 and SG-59), La Campana (LC-1 and LC-2) and Nahuelbuta (NA−1 and NA-2). DM, detergent mix added.